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An algorithm for numerical computation of the barrier transparency for the potentials 
surrounding Schwarzschild black holes is described for massless spinor particles. It is then 
applied to calculate the all configurations (including the contributions of twisted field 
configurations connected with Dirac monopoles) luminosity for the Hawking radiation 
from a Schwarzschild black hole. It is found that the contribution due to monopoles can 
be of order 22 % of the all configurations luminosity. 

1. Introductory Remarks 

While studying quantum geometry of fields on black holes the nontrivial topological 
properties of the latters may play essential role. The black holes can actually carry 
the whole spectrum of topologically inequivalent configurations (TICs) for miscella- 
neous fields, in the first turn, complex scalar and spinor ones. The mentioned TICs 
can markedly modify the Hawking radiation from black holes. Physically, the exis- 
tence of TICs should be obliged to the natural presence of magnetic U(N)-monopoles 
(with N > 1) on black holes though the total (internal) magnetic charge (abelian or 
nonabelian) of black hole remains equal to zero. Up to now, however, only influence 
of the TICs of complex scalar field on Hawking radiation has been studied more or 
less (see Ref. 1 and references therein for more details). The description of TICs for 
spinors was obtained in Ref. 2 but the detailed analysis of the TICs contribution to 
Hawking radiation requires knowledge of the conforming S-matrices which regulate 
the spinor particle passing through the potential barrier surrounding black hole. 
Those S'-matrices for the Schwarzschild (SW) black holes have only recently been 
explored in Ref. 3 which allow us to obtain an algorithm to calculate the S-matrix 
elements numerically since for physical results to be obtained one needs to apply 
the numerical methods. In the massless case the S'-matrices discussed are simpler 
to treat and the present paper will contain a description of the algorithm and will 
apply it to calculate the all configurations luminosity for massless spinor particles 
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for a SW black hole. The results obtained can serve as an estimate, in the first 
turn, of the electron-positron Hawking radiation and also for the neutrino one from 
the SW black holes. The more exact computation should take into account the 
particle masses and will require a more complicated algorithm for calculating the 
corresponding S'-matrices. The case of complex scalar field shows 1 , however, that 
the massless limit can be a good approximation to the more realistic massive case, 
in particular, when the particles masses are small enough. As a result, the task 
under consideration here makes sense. 

We write down the black hole metric under discussion (using the ordinary set 
of local coordinates t, r, i?, ip) in the form 

ds 2 = adt 2 - a^dr 2 - r 2 (d$ 2 + sin 2 ddip 2 ) ( 1 ) 

with a = 1 — 2M/r and M is the black hole mass. 

Throughout the paper we employ the system of units with h = c = G = 1, 
unless explicitly stated otherwise. Finally, we shall denote L 2 (F) the set of the 
modulo square integrable complex functions on any manifold F furnished with an 
integration measure while L^iF) will be the n-fold direct product of L,2(F) endowed 
with the obvious scalar product. 

2. Description of algorithm 

As was disscussed in Ref. 2, TICs of a spinor field on black holes are conditioned 
by the availability of a countable number of the twisted spinor bundles over the 
R 2 x S 2 -topology underlying the 4D black hole physics. From a physical point 
of view the appearance of spinor twisted configurations is linked with the natural 
presence of Dirac monopoles that play the role of connections in the complex line 
bundles corresponding to the twisted spinor bundles. Under the circumstances each 
TIC corresponds to sections of the corresponding spinor bundle E, which can be 
characterized by its Chern number neZ (the set of integers). Using the fact that all 
the mentioned bundles can be trivilized over the chart of local coordinates (t, r, ip) 
covering almost the whole manifold M 2 x § 2 one can obtain a suitable Dirac equation 
on the given chart for TIC \& with mass ito and Chern number n € Z that looks as 
follows 

Z> n * = iio*, (2) 

with the twisted Dirac operator T> n = i7 M V™ and we can call (standard) spinors 
corresponding to n — untwisted while the rest of the spinors with n^O should be 
referred to as twisted. Referring for details and for explicit form of V n to Ref. 2, it 
should be noted here that in L|(K 2 xS 2 ) there is a basis from the solutions of (2) 
in the form 

^ _ 1 c ia,« r -i ( ^i(r,w,A)<I> Am \ 
Am ~V2^ {F 2 (r,w,X)a^ Xm J ' [6) 

where <j\ is the Pauli matrix, the 2D spinor $>\ m = $^ TO (i?, <p) = ('I'lAm, ( E > 2Am) is 
the eigenspinor of the twisted euclidean Dirac operator with Chern number n on 
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the unit sphere with the eigenvalue A = ±y/ (I + l) 2 — n 2 while — I < m < I + 1, 
I > \n\. As said above, in the given paper we consider fi = and then the functions 
Fi t 2 obey the system of equations 



V^d r F 1 + + Fx = -icF 2 , 

y/M r F 2 + l^-$)F 2 = -icF 1 



(4) 



with c = u/yfa and a of (1). The explicit form of the 2D spinor $j m is inessential 
in the given paper and can be found in Ref. 2. One can only notice here that they 
can be subject to the normalization condition at n fixed 

(|<I'L.\ ro | 2 + |<I»o A „,|-)MniWiWv = I 
o o 

and these spinors form an orthonormal basis in L^S 2 ) at any n£Z. 

By passing on to the Regge- Wheeler variable = r + 2Mln(r/2M — 1) and 
by going to the quantities x = r*/M,y = r/M,k = loM, we shall have x = y + 
2 ln(0.5y — 1), so that y(x) is given implicitly by the latter relation (i.e., — oo < x < 
co, 2 < y < oo) with 

y' = dy/dx =l-2/y=(y- 2)/y = a (x) (5) 

and the system (4) can be rewriten as follows 

E[ + a\E\ = b\E 2 , , . 

E' 2 + a 2 E 2 = b 2 E x [0) 

with^i =Ei(x,k,\)=F + (Mx),F+(r*) = F 1 [r(r*)],E 2 = E 2 (x,k,\) = iF_(Mx),F_ 
F 2 [r(r*)] and 

aia =2~f ± y^' (?) 
h, 2 = Tk. (8) 

To evaluate luminosity of the Hawking radiation for spinor particles it is necessary 
to know the asymptotics of the functions E\, 2 at x — > +oo. As was shown in Ref. 
3, the latter asymptotics look as follows 

Akx 



E 1 - V-fcsn (k, \)e lkx x -> +oo , (9) 

^ MMM) ^ ,^ +00) (10) 

where Sn(fc, A) is an element of the S'-matrix connected with some scattering prob- 
lem for the Schrodinger-like equation 

u" + k 2 u = qu (11) 
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with potential 

y{x) 



q{x,X) = -TTTV / ao 
y z (x) 



(12) 



Xy(x) 

In its turn, the correct statement of the mentioned scattering problem for Eq. (11) 
consists in searching for two solutions u + (x,k,X), u~(x,k,X) of the equation (11) 
obeying the following conditions 



u + (x, k, A) 



e tkx + s 12 (k,X)e- tkx +o(l), x^-oo, 
su(k, X)e lkx + o(l), x — > +oo, 



u-{x,k^) = { s ^ e ~ ik ;+Sl\ m * ^ 7°°' (13) 
(e tKx + s 2 i(k, X)e lhx + o(l), a; — > +oo. 

It is seen that there arises some 5- matrix with elements Sij,i,j = 1,2. After this 
one can obtain the luminosity L(n) with respect to the Hawking radiation for TIC 
with the Chern number n in the form (for more details see Ref. 2 and 3) 



L(n)=A^f^2(l + l) [ 



±A l=\n\ 



, £ 2(, + Wy^>'V (U, 

H™l 

where .4 = (^) 1/2 ~ 0.251455 • 10 55 erg • s" 1 • M" 1 and M in g. Luminosity 
L(n) can be interpreted, as usual, 1,2 as an additional contribution to the Hawking 
radiation due to the additional spinor particles leaving the black hole because of the 
interaction with monopoles and the conforming radiation can be called the monopole 
Hawking radiation. 2 

Under this situation, for the all configurations luminosity L of black hole with 
respect to the Hawking radiation concerning the spinor field to be obtained, one 
should sum up over all n, i. e. 

oo 

L = J2 L(n) = L(Q) + 2j2 L(n) (15) 

raSZ n=l 

since L{— n) = L(n). 

It should be emphasized that in (14), generally speaking, sn(/c, A) ^ sn(fc, —A) 
and that is clear from the form of potential q of (12) which satisfies the condition 3 

+oo 

J \q(x, X)\dx < oo . 



As an example, in Figs. 1, 2 the numerically computed potentials q(x, A) with A > 
and A < are shown. 
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Fig. 1. Typical potential barrier for A > 0. 



Fig. 2. Typical potential barrier for A < 0. 
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One can notice that when A < 0, q(x, A) > at any x € [— oo, oo] while at A > 0, 
q(x, A) can change the sign. 

Obviously, for neutrino there exists only L(0) since neutrino does not interact 
with monopoles (when neglecting a possible insignificantly small magnetic moment 
of neutrino). Also it is evident that for numerical computation of all the above 
luminosities one needs to have some algorithm for calculating Sn(fc, A). This can 
be extracted from the results of Rcf. 3. Namely 

s n (k,X) =2ik/[f-(x,k,X),f+(x,k,X)} , (16) 

where [,] signifies the Wronskian of functions f~,f + , the so-called Jost type so- 
lutions of Eq. (11). In their turn, these functions and their derivatives obey the 
certain integral equations. Since the Wronskian does not depend on x one can take 
the following form of the mentioned integral equations 

f-(x ,k,X) = e- ikx ° + ^ J sm[k(x -t)}q(t,X)f-(t,k,X)dt, (17) 

— oo 
x a 

{f-)' x {x ,k,X) = -ike- ikxo + J cos[k{x -t)]q(t,X)f-(t,k,X)dt, (18) 

— OO 

+oo 

f + (x ,k,X) =e lkx ° +1 J sm[k(x -t)}q(t,X)f+{t,k,X)dt, (19) 

Xq 

+oo 

(f + y x (x ,k,X) =ike lkxa + J cos[k{x -t)]q{t,X)f + {t,k,X)dt. (20) 

The point xq should be chosen from the considerations of the computational con- 
venience. The relations (17)-(20) can be employed for numerical calculation of 
su(k, A). 

Finally, we should touch upon the convergence of the series (15) over n. For this 
aim we denote 



e 87Tk _|_ I 





yit-Lij/^n OO 

i (±A) = ca(±A) + c /2 (±A)= J + J , (21) 



and we represent the coefficients c/(±A) in the form (omitting the integrand) 

(lni)/2ir 

r/( ' 

(lni)/2ir 

so that L(n) of (14) is equal to L\{n) + L,2(n) respectively. Also it should be noted, 
as follows from the results of Ref. 3, we have for barrier transparency T(k, A) = 
|sn(fe, A| 2 that T(0, A) = and at k -> +oo 

| Sll (fc,A)| 2 =r(fc,A) = l + 0(fc- 1 ). (22) 
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Further we have the asymptotic behaviour for large I and k « I 



-Vn+l+l 

\ Sll (k,l,n)\=C == fe[l + o(l)] (23) 
Vn + 1 + 1 

with some constant C. To obtain it we should work using the equations (17)-(20) 
by conventional methods of mathematical analysis so that the detailed derivation 
of (23) lies somewhat out of the scope of the present paper and will be considered 
elesewhere. Using (21) and (23) we find 

e -2V«+i+l /-(ln/)/2ir ^3 

c/i(±A) < B ; — - / -z-r — -dk 

n ; - n + l + 1 J e^ k + l 

e -2Vn+l+l roc .3 £ -2Vn+l+l Q\/-(4) 



n + l + 1 J e^ k - 1 n + l + 1 (8tt) 4 

with some constant B, where we used the formula (for natural p > 0) 4 



00 

J i"dl 




£^-=p!C(P + l) 



with the Riemann zeta function ((s), while C(4) = 7r 4 /90. As a consequence, one 
may consider 

OO . oo 

/p-2Vn+(+l f / T)\ 

2 ( / + i ) - TTTrrf /~2 y e - 2t (t- 7 )^ 

(24) 



where we employed the asymptotical behaviour of the incomplete gamma function 4 

oo 

T(a, x) = y e-H a - x dt ~ a^-V*, a; -» +oo , 

so that the series ^i( n ) is evidently convergent. 
At the same time due to (22) 

CO 

c l2 (±\)<D J ke~ 8wk dk ~ (25) 

(ln/)/27r 

with some constant D. This entails 



oc 



, . f 2(f + l)lnZ „ Inn 
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and the series J2T ^(n) is also convergent. Under this situation we obtain that 
the series of (15) is convergent, i. e., L < +00. Consequently, we can say that all 
the further computed luminosities are well defined and exist. 

3. Numerical results 

In view of (22), we have 

00 00 

r mxjkdk r kdk 1 „ , 1 (97 , 

J e 8 - fc + l 7 e 8wfe -l (8tt) 2 U j 6 • 8 2 ' { ' 



whilst the latter integral can be accurately evaluated using the trapezium formula on 
the [0, 3] interval. Having confined the range of k to the [0, 3] interval, accordingly, 
it is actually enough to restrict oneself to < I, n < 15-20 when computing L(n), L. 
We took < I < 20, < n < 15. Since the Wronskian of (16) does not depend 
on x, the latter should be chosen in the region where the potentials q(x, ±A) arc 
already small enough. Besides, potentials q are really equal to when x < -30. We 
computed T(k, A) according to (16) at x = x = 300, where were obtained from 
the Volterra integral equations, respectively, (17) and (19). For this aim, according 
to the methods developed for numerical solutions of integral equations (see e. g., 
Ref. 5), those of (17) and (19) have been replaced by systems of linear algebraic 
equations which can be gained when calculating the conforming integrals by the 
trapezium formula, respectively, for the [-30, 300] and [300, 400] intervals. The 
sought values of were obtained as the solutions to the above linear systems. 
After this the derivatives (Z 1 * 1 )^ were evaluated in accordance with (18) and (20) 
while employing the values obtained for f ± . The typical behaviour of T(fc, A) is 
presented in Fig. 3. 



Fig. 3. Typical behaviour of barrier transparencies. 



Finally, Fig. 4 represents the untwisted L(0)/A and the all configurations L/A 
luminosities with A of (14) as functions of k. The areas under the curves give the 
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corresponding values of L(0)/A and L/A 

L(0)/A = 0.514490 • , 

L/A = 0.658571 • 10~ 3 , 
L(0)/L = 0.781222, 

so that contribution owing to Dirac monopoles amounts to 21.8778% ~ 22% of L. 
For inquiry, the similar contribution for scalar particles was of order (17-18) %. 1 



Fig. 4. Untwisted and all configurations luminosities. 



4. Concluding remarks 

As was mentioned above, the results obtained can serve as an estimate, in the 
first turn, of the electron-positron Hawking radiation and also of that for neutrino 
from the SW black holes. Clearly, for neutrino there exists only L(0) since neutrino 
does not interact with monopoles (when neglecting a possible insignificantly small 
magnetic moment of neutrino). The corresponding computation for other funda- 
mental fermions (/i ± -mesons and r ± -leptons) and also the more exact computation 
for e^- -particles should take into account the particle masses and will require a more 
complicated algorithm for calculating the corresponding S'-matrices. We hope to 
get the numerical results in the massive case elsewhere. 
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